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Abstract 

In this paper the local order of convergence used in iterative methods to solve nonlinear systems of 
equations is revisited, where shorter alternative analytic proofs of the order based on developments of 
multilineal functions are shown. Most important, an adaptive multi-precision arithmetics is used hereof, 
where in each step the length of the mantissa is defined independently of the knowledge of the root. 
Furthermore, generalizations of the one dimensional case to m-dimensions of three approximations of 
computational order of convergence are defined. Examples illustrating the previous results are given. 
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1 Introduction 

Well-known analytic techniques [E]-[8] to prove the local order of convergence of iterative methods to solve 
single nonlinear equations are generalized to systems of nonlinear equations. More precisely, a new proof of 
the local order for an iterative method to solve the system of equations F(x) = 0, where F : W n — > R m , is 
presented. The basic tools used are the formal developments of the function F, its inverse and its derivatives 
in power series. The vectorial expression of the error equation obtained carrying out this procedure, as we 
will see later on, is 

e n+1 = H (F'(a), F"(a), ...)<+ O , 
where a is a simple root of F(x) = 0. 

The preceding technique to prove the local convergence is illustrated with several examples in which 
generalizations of the unidimensional case to m-dimensions are performed using the following definitions 
of computational order of convergence: 

• Computational order of convergence (COC) in Weerakoon et al. (2000, [T9"]V 

• Approximated order of convergence (ACOC) in Hueso et al. (2009, |12j). 

• Extrapolated order of convergence (ECOC) in Grau et al. (2009, [9]). 
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2 Notation and basic results 



To obtain the vectorial equation of the error we need some known results that for ease reference are included 
in the following. Let F : D C R m — > K m be sufficiently differentiable (Frechet-differentiable) in D, and 
therefore with its differentials continuous. If we consider the fcth derivative of F at a € R m , we have the 
/c-linear function 



k 

•■x ■■■ xl 



(hi,...,h k ) .— > F^(a)(h u ...,h k ). 
That is, F^ k) (a) (hi, . . . , h k ) e R m . It has the following properties: 
PI. F^(a) {hi,..., h k -i, • ) G if (M m , M m ) = if (R m ). 

P2. F^(b) (/icr(i), ■ • ■ , ^ CT (/s)) — F^ k \a) (hi, ... , ftfe), where cr is any permutation of the set {1, 2, . . . fc}. 
Notice that from PI and P2 we can use the following notation: 

Nl. f!"(«)C'i M = f ( ' ) («)*i"^- 

N2. F<*> (a) h k - x (a) ^ = F<*> (a) (a) /i fe +^" 1 . 
Hence, we also can express F^ k '(d) (hi, . . . , h k ) as 

F {k \a) (hi,...,h k -i)h k =F^\a) {h u ...,h k -i) h k -ih k = ... = F (k \a) hi ■ ■ ■ h k . 

On the other hand, for any q = a + h 6 M m lying in a neighborhood of a simple zero, a G R m , of the system 
F(x) = we can apply Taylor's formulae and assuming that there exists T = [F'(a)]~ , we have 

3 

h + ^ A k h k + 4 



where 



F(a +h) = F'(a) 
1 



k=2 



(1) 



A k = — T F (K, (a), k>2, and 4 = 0(h*). 



Note that since F^ (a) G if 
express the differential of first order as: 



fc 

x xl 



l m Tip m 



and T G if (M m ), then A fe ft, fe G R m . Moreover, we can 



/ + ^ M fe Z^- 1 + 3 



k=2 



F'(a + h) =F'(a) 

k 

where I is the identity and M k = k A k . Therefore, M k h^ 1 G if (K m ). From {2}, we get 

[f'(q + ft)] -1 = [/ - g 2 ft + g 3 ft 2 + o 3 ] r, 

f Q 2 = 



(2) 



(3) 



where 

Hereof we have used the following notation: M k i h k+e ~ 2 = M k h k ~ x Mg h 1 " 1 , if k ^ I, and Mf instead of 
M«, ]£k=£. 



Q 3 = Ml - M 3 
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In general, if [F'id)}^ 1 = [I - Q 2 S + Q 3 6 2 ] T and F(t) = F'(a) [r + A 2 r 2 + A 3 r 3 ] , with the same 
notation, we obtain: 



[F'(d)} 1 F(t) = t+(A 2 t 2 -Q 2 5t) + (A3 t 3 - Q 2 A 2 St 2 + Q 3 S 2 t) 



(4) 



We close this section applying the preceding to the Newton's method. Let x € K m , Newton's method is 
given by 

z = x- \F'{x)Y x F(x). (5) 

Putting q = x in (JTJ) and ([3]) we obtain F(x) and [F 1 in powers oi e = x — a. The expression of 

= z — a in terms of e is build up subtracting a to both sides of ([5]) . Namely, 



E = e - [F\x)]-' F{x) = T 2 e 2 + T 3 e A + (9 4 



(6) 



where T 2 = M 2 /2 = A 2 , and T 3 = (4M 3 - 3Mf)/6 = 2 (A3 - These values agree with the classical 
asymptotical constant in the one dimensional case. 

Without to use norms we can define the local order of convergence for one-step iterative method as follows. 

} ;,; !!; 1 such 



The local order of convergence is p € N if there exists a p-linear function K 6 Jz? yBJ n x — ■ x I 
that 

E = K eP + O p+1 , (7) 

p 

where e p is (e, e) and O p +i = 0(e p+1 ). 



3 Variants of Newton's method 

In the following, using the results presented in the previous section, three known variants of Newton's method 
with local order of convergence equal three, are analyzed. We explicitly give their vectorial error equation in 
which it appears a 3-linear application instead of the asymptotical error constant used in the one dimensional 
case. 



3.1 Arithmetic mean Newton's method 

The first variant ([Hj,[19]), that substitutes the derivative of F(x) by the arithmetic mean of the derivatives 
of F at the points x and z is given by 

X = x- 2[F'(x)+F'(z)}- 1 F{x). (8) 
From (J6j) and the developments of F'(x), F'(z), we derive the development of [F'(x) + F'(z)] , Indeed, 



F'{x) 
F'(z) 

where T 2 = M 2 /2. From (J3j) we get 
[F'(x)+F'(z)}- 1 



F'(a) [I + M 2 e + M 3 e 2 ] + 3 , 
F'{a) [I + M 2 E] + 0(E 2 ) 
F'(a) [l + M 2 T 2 e 2 ] +0 3 , 

[F'(a) {21 + M 2 e + (M 3 + M 2 T 2 ) e 2 }] ^ 



I- (M 2 e + M 3 e 2 ) 



r +o 3 



(9) 
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Taking into account that = k Ak, fc = 2, 3 , subtracting a from ©, and applying Q we have 

X-a = {-A 3 ' " 1 ' 1 ' 1 ' 1 ■ - ; 
T 



^M 2 A 2 + ^M 3 j - ' (>, 



o 4 



(10) 



Note that we get the same expression as the one obtained for the one dimensional case in [18], [19] and [6]. 
3.2 Harmonic mean Newton's method 

The second variant ( [18] , [IT] , [15] ) , that substitutes the derivative of F(x) by the harmonic mean of the 
derivatives of F at the points x and z is given by 



X= x- - [F'ixy'+F'iz)- 1 } F(x). 

From the developments of F'(x), F'(z) we derive the developments of F'fa)^ 1 and F'(z) -1 . Namely, 

F' (x)- 1 = [I- M 2 e + (M| - M 3 ) e 2 ] r + 3 , 

F'iz)- 1 = [I-M 2 E]T + 0{e 4 ) 
1 



Mn e 2 



r + o 3 . 



From (IT2l we have 



Finally, subtracting a from (fTTj) .taking into account (IT3"|) and applying yields 



1 /l 



-A 3 + -M 2 A 2 -- 



-A3 e 3 + 4 . 



Mo - Mp 



o 4 



This result agrees with the ones obtained in [18] . [TTj and [15] . 



(11) 



(12) 
(13) 



(14) 



3.3 Frozen derivatives in Newton's method 

The third iterative method presented in this section is derived independently in the works of Shamanskii 
([17], 1967) and Potra et al. ([16], 1989). It is defined by 

X = z- F'ix^Fiz), (15) 

where z = x— F 1 (x)^ 1 F(x) is a Newton point. This method is a modification of Newton's iterative function 
where F is computed in the second step without evaluating F'. So, we can consider that the derivative is 
frozen in this second step ([H],[T4]-[2]). Note that only one computation of the inverse function is needed 
and it is only necessary one LU decomposition. 

Subtracting a from (|15p and applying Q we obtain the following vectorial error equation: 

X-a = E-F'(x)- 1 F(z) 

= T 2 e 2 + T 3 e 3 -[I- M 2 e + 2 ] [A\ e 2 + 2 (A 3 - A 2 2 ) e 3 + 4 ] + O4 

= 2A 2 > 3 + 4 - (16) 



4 



4 Main numerical results 



In this section we introduce two variants of the definition of COC which are independent of the knowledge of 
the root and we give some ways to compute with adaptive multiprecision arithmetics and a stopping criteria. 
Moreover, we apply all the results obtained yet to several examples. 



4.1 Theoretical concepts 

Let x„_2, x n -i, x n and be the last four consecutive iterations of the sequence {x n } n >o, where x n £ R. 
The definitions of Computational Order of Convergence (COC), p n [19] . Approximated Computational Order 
of Convergence (ACOC), p n [12], and Extrapolated Computational Order of Convergence (ECOC), p n [9], 
are the following 

Definition 1 The COC, ACOC and ECOC of a sequence {x n } n >o is defined by 

_ _ ln|e» + i/e„| „ _ In | e n+ i/e„| „ _ In | e» + i/e»| 

In |e„/e„_i| ' ™ ln|e„/e„_i|' " In | e„/e„_i| 

respectively, where e n — x n — a, e n — x n — x„_i , e n — x n ~ a n and a n = x n — ^'i^ 1 ^ , ri> 2, where A 
is the forward difference operator, Axk — Xk+i — x^. 



One of the main drawback of the COC is that it involves the exact root a, which in a real situation it is 
not known a priori. To avoid this, we have introduced these two variants of COC, that do not use the exact 
root. 

In numerical problems where a huge number of significant digits of the solution is needed it is required the 
use of methods with a high order of convergence together with an adaptive arithmetics that is to update the 
length of the mantissa at each step by means of the formula 

Digits:= [px (- log | e„| + j)] , (17) 

where p is the order of convergence of the method and [x] denotes the integer part of x. Notice that the 
length of the mantissa is increased approximately by the order of convergence p. We have numerically 
checked the value of j, by varying it between 1 and 5, in order to have enough accuracy in the computation 
of the iterates {x n } n >o- We have realized that the minimum value that guarantees all the significant digits 
required is, in almost cases, j — 2. Consequently, hereof we consider j — 2 in formula (ITTl) except if another 
value is explicitly given. In addition, to compute e„, e n or e„ with an appropriate number of figures, using 
Definition 1, we must to enlarge the mantissa in the computation of x n+ i,x n ,x n _i, . . . with at least four 
additional significant digits. 

Two relationships, one between e„ and e„, and the other between e n and e„ are given in |10j . Namely, 

e „« C^h (J^-Y /iP 1] and e„ « C ~ef'^' l \ (18) 
\e n -ij 

Notice that for updating the adaptive arithmetic process (fTT]) it is necessary to know the exact root a. In 
this case the following stopping criteria is applied: 

|e„| = \x n - a\ < 0.5 • 10"", (19) 
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where r\ is the number of correct decimals and 0.5 • 10~'' is the required accuracy. The result given in (|18[) 
allows us to substitute the error in (1171) by an expression that does not involve the exact root. Indeed, we 
implement the following adaptive multi-precision arithmetic schemes: 



Digits 

where S n — 



P-1 



x (-log|<5 n | + j) 



P" 

2/5-1 



x (-log|e n | + j) 



(20) 



e ra -i 



or Digits 

Moreover, from (fT8|) we propose the following stopping criteria, instead of p^|) : 

5 n < 0.5- IQ-'^p-V/p 2 or |e„| < 0.5 • lO^ (2 "~ 1)/p2 . (21) 



4.2 Examples 



We generalize the preceding definition and techniques for solving seven systems of nonlinear equations using 
the Maple computer algebra system. We use the norm || • ||oo instead of the absolute value in ([!?]) - ([2"TT) . 

We have computed the solution of each system for the same set of initial approximations x l , i = 1,2,3, 
which have been chosen close to the root a using the eucliden distance di = \\x l — ct 1 1 2 , and taking into 
account the value of Di = ||-F(a:o)lloo- 

Depending on the computational order of convergence used, COC, ACOC or ECOC, the iterative method 
was stopped when the condition (|19l) or (|2ip is fulfilled. Note that in all cases 77 = 2800 and we also obtain 
||F(xfe)|| < 0.5-10~ ,? . Finally, we choose j, see (|T7|) and ([20]) . such that ||efc||, ||efe|| or ||e&|| respectively, and 
||F(xfe)|| have three significant digits al least. 

Tables 1-7 show, for each method and each function, the number, k, of iterations needed to compute the root 
to the level of precision described. Note that independently of using (flTf or (|20[) the number of necessary 
iterations to get the desired precision is the same. In addition, in the sixth, eighth and tenth column, it 
is shown an error bound for the corresponding Computational Orders of Convergence (COC, ACOC and 
ECOC), given respectively by p k -i = P ± Ap fe _i, p k = p± Ap k , and p k — p± Ap k . 

4.2.1 Example 1 

We begin with the system iq (x) = defined by 

f e x - 2 = 0, 
[ sin(2y - x) = 0. 

The roots of -Fi(a;) = are (x, y) = (In 2, kir + In \/2)* with k € Z. We study the convergence of iterative 
methods previously presented towards the root a = (In 2, lnV^)* (0.6931471806, 0.3465735903)' . 
The initial approximations of the four methods are x\ = (1, 0)* with d\ = 0.347 and D\ — 0.841; 
x\ = (0.6, 0.3)' with d 2 = 0.0931 and D 2 = 0.178; and x% = (0.7, 0.35)* with d 3 = 0.00685 and 
D 3 = 0.0137. 

Note that in the results shown in Table [T] in the case marked with (*) it is necessary to take j — 4 for x\ 
and Xq in (fl"7|) and ((20|) because for j < 4 we loose some significant figures in the computation of iq (x k ) . In 
the case marked with (*), for the same reason, we take j = 8 for Xq. 
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Tabic 1: Numerical results for system Fx(x) = . 



















coc 










ECOC 










ACOC 








Method 


Xo 


k 




Xk) 


'IL 


lie* 


-l 


II oo 


Ap 


fc-i 




II 


e*|l„ 




Ap fe 














NM 


Xq 


12 


1.02 • 


10" 


-3429 


3.19 ■ 


10 


-1715 


4.03 • 


10" 


-15 


1.19 


■ 10" 


2572 


3.78 ■ 


10" 


-19 


3.99 


■ 10" 


-858 


6.74 


• 10 


-20 


p = 2 


T 2 


12 


2.95 • 


10" 


-5428 


1.72 ■ 


10 


-2714 


5.17 • 


10" 


-16 


1.59 


■ 10" 


4071 


4.16 ■ 


10" 


-20 


9.27 ■ 


10" 


1358 


1.44 


• 10 


-25 


(*) 


x 3 


11 


8.63 • 


lo- 


-5050 


2.94 ■ 


10 


-2525 


1.91 • 


10" 


-15 


1.13 


■ 10" 


3787 


7.61 • 


10" 


-19 


3.83 ■ 


10" 


1263 


3.32 


• 10 


-25 


AMN 




9 


1.04 • 


ur 


-7959 


4.93 ■ 


10 


-2832 


2.69 • 


10" 


-17 


1.91 


■ 10" 


4422 


1.40 ■ 


10" 


-13 


1.65 ■ 


10" 


1769 


2.34 


• 10 


-35 


p = 3 


T 2 


8 


3.46 • 


10" 


-8274 


1.73 ■ 


10 


-2758 


9.86 • 


10" 


-18 


3.73 


■ 10" 


4597 


1.83 ■ 


10" 


-17 


2.15 ■ 


10" 


1839 


9.74 


• 10 


-36 




t 3 

'''0 


7 


3.13 • 


10" 


-5256 


1.68 ■ 


10 


-1752 


9.86 • 


10" 


-18 


1.64 


■ 10" 


2920 


4.73 ■ 


10" 


-18 


9.78 ■ 


10" 


1169 


2.01 


• 10 


-35 


HMN 


Xq 


9 


6.80 • 


10" 


-7162 


1.01 ■ 


10 


-2387 


5.80 • 


10" 


-20 


3.25 


■ 10" 


3979 


2.05 


10" 


-19 


3.23 ■ 


10" 


1592 


5.00 


• 10 


-35 


p = 3 


T 2 

'''0 


7 


2.83 • 


10" 


-3451 


1.19 ■ 


10 


-1150 


3.26 • 


10" 


-20 


1.26 


■ 10" 


1917 


2.63 


10" 


-21 


1.06 


■ 10" 


-767 


5.35 


• 10 


-35 




r 3 


7 


7.52 • 


10" 


-5912 


7.67 ■ 


10 


-1971 


5.00 • 


10" 


-18 


1.30 


■ 10" 


3284 


3.06 


10" 


-21 


1.70 ■ 


10" 


1314 


1.72 


• 10 


-31 


FDN 


x 1 


8 


2.52 • 


10" 


-4653 


1.36 ■ 


10 


-1551 


3.12 • 


10" 


-32 


2.98 


■ 10" 


-862 


5.02 ■ 


10" 


-16 


9.75 ■ 


10" 


1035 


2.52 ■ 


10" 


-117 


p = 3 


X 2 




8 


5.46 • 


10" 


-7653 


1.76 ■ 


10 


-2551 


2.42 • 


10" 


-17 


4.39 


■ 10" 


4252 


3.10 


10" 


-18 


2.49 ■ 


10" 


1701 


9.24 


• 10 


-36 


0) 


r 3 


7 


9.43 • 


10" 


-5065 


9.81 ■ 


10 


-1689 


9.63 • 


10" 


-18 


3.57 


■ 10" 


2814 


1.35 ■ 


10" 


-17 


3.64 ■ 


10" 


1126 


4.06 


• 10 


-35 



4.2.2 Example 2 

The second example involves two quadratic polynomials. Namely, F 2 (x) = 0, defined by 

J x 2 - 4a; + y 2 = 0, 
| 2x + y 2 -2 = 0. 



Its solutions are (x,y) = - y/7, ±yj - 4 + 2^/7 j and (x, y) = ^3 + ■yff, ±yj - 4 - 2^7 J . We test the 

convergence of the methods towards the root a = ^3 - V7, \J — 4 + 2 V^^) ~ (0.3542486889, 1.136442969)* . 

The initial values are x\ = (-1, 0.4)' with d x = 1.354 and D x = 5.16; x 2 = (0, 1)' with d 2 = 0.354 and 
D 2 = 1.0 ; and = (0.3, 1.1)* with d 3 = 0.0542 and D 3 = 0.19 . 

Note that in Table [2] the method HMN is of 4th order according to (IT6|) . Setting in (*) p = 4 the results 
obtained for ECOC and ACOC were excellent for the three initial conditions with j = 2, but for COC it 
was necessary to take j — 20. Hence in the computations shown in Table[2]was used j = 20. In (*), for NM 
method j = 3 was used. 

4.2.3 Example 3 

A system of equations involving cubic polynomials, F 3 (x) — 0, namely 

j x 3 - 3xy 2 -1 = 0, 
\ 3x 2 y-y 3 + 1 = 0, 

is analyzed. 
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Table 2: Numerical results for system F 2 (x) = . 



















coc 










ECOC 










ACOC 








Method 


Xo 


A- 




Xk) 


'Hoc 


IK 


-l 


Hoc 


Ap 


fc-i 




11 


e*ll„ 




Ap fe 










Ap fe 




NM 


Xq 


14 


4.33 • 


10" 


-3666 


2.08 ■ 


10 


-1833 


1.68 • 


10" 


-14 


3.78 


• 10" 


1375 


2.35 ■ 


10" 


-15 


3.03 


■ 10" 


-917 


4.98 ■ 


10" 


-40 


p = 2 


T 2 


12 


3.34 • 


10" 


-3778 


1.83 ■ 


10 


-1889 


2.38 • 


10" 


-14 


3.43 


• 10" 


1417 


2.98 ■ 


10" 


-15 


5.33 


• 10" 


-473 


3.26 


10" 


-40 


(*) 


a: 3 

Xq 


11 


3.46 • 


10" 


-3375 


5.88 ■ 


10 


-1688 


1.88 • 


10" 


-13 


9.46 


• 10" 


2532 


5.02 ■ 


10" 


-12 


1.27 


■ 10" 


-422 


3.16 


10" 


-32 


AMN 




9 


6.91 • 


10" 


-5380 


1.16 ■ 


10 


-1793 


8.54 • 


lo- 


-11 


3.45 


• 10" 


2989 


3.96 


10" 


-18 


2.97 ■ 


10" 


1196 


6.45 ■ 


10" 


-36 


p = 3 


t 2 


8 


2.40 ■ 


10" 


-7591 


8.17 ■ 


10 


-2531 


8.68 ■ 


ur 


-13 


8.90 


• 10" 


4218 


4.23 


10" 


-11 


1.09 ■ 


10" 


1687 


1.40 ■ 


10" 


-25 




Xq 


7 


3.95 ■ 


10" 


-3817 


9.64 ■ 


10 


-1273 


2.77 ■ 


10" 


-13 


5.45 


• 10" 


2121 


7.54 ■ 


10" 


-17 


5.65 


■ 10" 


-849 


5.59 ■ 


10" 


-35 


HMN 


x q 


7 


4.33 • 


10" 


-3666 


5.20 ■ 


10 


-1230 


1.06 • 


lo- 


-16 


2.81 


• 10" 


1604 


2.76 


10" 


-17 


4.08 


■ 10" 


-688 


1.20 ■ 


10" 


-28 


p = 4 


T 2 
Xq 


6 


3.34 • 


10" 


-4067 


6.45 ■ 


10 


-945 


1.45 • 


ur 


-38 


2.51 


• 10" 


1653 


6.29 


10" 


-12 


3.89 


■ 10" 


-709 


1.58 ■ 


10" 


-28 


(*) 


r 3 


6 


2.32 • 


10" 


-6750 


5.88 ■ 


10 


-1688 


3.27 • 


10" 


-47 


1.20 


• 10" 


2953 


2.81 ■ 


10" 


-14 


2.04 ■ 


10" 


1266 


3.27 ■ 


10" 


-47 


FDN 


■'-o 


10 


1.83 • 


10" 


-5372 


2.75 ■ 


10 


-1791 


8.88 • 


lo- 


-13 


3.94 


• 10" 


2985 


1.29 ■ 


10" 


-19 


1.43 ■ 


10" 


1194 


1.72 ■ 


10" 


-34 


p = 3 


X 2 


8 


9.36 ■ 


10" 


-6199 


1.02 ■ 


10 


-2066 


3.78 ■ 


ur 


-14 


3.50 


• 10" 


3444 


6.42 ■ 


10" 


-20 


3.43 ■ 


10" 


1378 


1.02 ■ 


10" 


-34 




r 3 
x 


7 


5.20 ■ 


10" 


-3475 


8.39 ■ 


10 


-1159 


2.28 • 


10" 


-13 


5.44 


• 10" 


1931 


2.94 ■ 


10" 


-17 


6.48 


• 10" 


-773 


1.15 ■ 


10" 


-34 



Its solutions are (x,y) = ^5 + 3^/3^ • (3' 2 = ( 5 + 3 V / 3) ' [~^ 1 '2I and 

(x,y) = (— 2~ 1 / 3 , — 2 -1 / 3 ). We test the convergence of the methods towards the root a = + 3\/3 

(4-$ 



(-0.2905145555, 1.084215081)* . 

The initial values are = (-1, 2)* with d x = 0.916 and D 1 = 10.0; x 2 = (-0.1,1.4)* with d 2 = 0.316 
and D 2 = 1.702; and xl = (-0.3,1.1)* with rf 3 = 0.0158 and D 3 = 0.062. Notice that in TableHall the 
iterative methods give the same values of Ap k _ 1 and Apk ■ 

4.2.4 Example 4 

We present an example consisting in the computation of the complex root of the equation e z = z with the 
smallest imaginary part. Here we solve the nonlinear system of equations 

{e x cosy = x, 
e x siny = y, 

that is obtained by considering the real and imaginary parts of the original equation e x+zy = x + iy. 

We apply the methods defined in previous sections setting = (0,2)* with d\ = 0.6627 and D\ = 1.091; 
xl = (0.2, 1.1)* with d 2 = 0.2372 and D 2 = 0.354 ; and x% = (0.3, 1.3)* with d 3 = 0.03723 and D 3 = 0.0611 . 

In Table|4] in the cases (*) and (*) was necessary to use j = 3 and j = 4 respectively, and again Ap k _ 1 and 
Apk agree. This example can be found [7]. 
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Table 3: Numerical results for system F 3 (x) = 



Method 



COC 



ECOC 



\\F 3 (^)\\ C 



[K-i|| c 



lie* ||, 



Ap fe 



ACOC 

Mlco 



NM 
p = 2 



14 9.85-10 
13 3.52-10 
11 3.97-10 



-3759 

-4507 
-3665 



5.43 ■ 10 -1880 

2.85- 10- 2244 
3.26 ■ 10" 1833 



9.00 • 10" 



3.68 • 10" 



6.29 • 10" 



2.82 • 10" 2819 
4.57 • 10" 3381 
1.77- 10" 2749 



2.60 ■ 10" 



5.59 ■ 10" 



2.08 ■ 10" 



2.10 
1.60 ■ 
5.46 



10" 940 
10 -1127 

10 -917 



9.00 ■ 10"° 
3.68 • 10" 4 
6.29 ■ 10~ 6 



AMN 
p = 3 



x 







9 3.20-10 
8 8.73 ■ 10 
7 2.23-10 



-4152 
-3422 
-3841 



9.37 ■ icr 1385 
2.54- 10" 1141 
3.60 ■ 10" 1281 



1.71 • 10" 
4.38 ■ 10" 
1.03 ■ 10" 



5.52- 10" 2307 
1.89- 10" 1901 
4.08 ■ 10" 2134 



3.88 ■ 10" 
5.50 ■ 10" 
3.78 ■ 10" 



1.97 
4.32 
2.25 



1Q -923 

10" 761 
■ 10" 554 



1.71 • 10" 
4.37 • 10" 
1.03 ■ 10" 



HMN 
p = 3 



x\ 9 5.18-10 
x 2 a 8 2.12 ■ 10 
xlo 7 1.64-10 



-7826 
-5383 
-4740 



4.95 ■ lO" 2609 
7.58 ■ 10" 1795 
1.20 ■ 10" 1580 



4.43 ■ 10" 
3.13 ■ 10" 
5.52 ■ 10" 



5.76 ■ lO" 4348 
8.65- 10" 2991 
3.78 ■ 10" 2633 



5.85 ■ 10" 
1.27 ■ 10" 
4.56 ■ 10" 



8.41 ■ 
4.55 ■ 
2.47 ■ 



10 -1740 

10" 1795 

1Q -1054 



4.43 ■ 10" 
3.13 • 10" 
5.52 • 10" 



FDN 
p = 3 



x, 



8 3.40 • 10 
8 2.38-10 
7 1.92 • 10 



-3218 
-2869 
-3591 



1.58 ■ 10" 1073 

2.87 ■ 10" 957 
5.43 ■ 10" 1198 



8.07 • 10" 4 
9.83 ■ 10" 4 
1.10 ■ 10" 4 



2.19 ■ 10" 178S 
8.18 ■ 10" 1B9B 
6.60- 10" 1996 



4.39 • 10~ 3 
2.18 ■ 10" 3 
2.76 ■ 10" 3 



8.23 
2.85 
8.45 



■ lO" 718 

■ 10" 638 

■ 10" 799 



8.07 • 10" 4 
9.83 • 10" 4 
1.10 • 10" 4 



4.2.5 Example 5 

A system F 5 (x) = involving three nonlinear equations denned by 

xyz = 1, 
x + y - z 2 = 0, 
x 2 + y 2 + z 2 = 9, 

is studied. We analyze the convergence of the methods towards the root 

a « (2.14025812200, -2.09029464225, -0.22352512107)*. 

The initial values are x\ = (1.0, -1.0, 0.1)' with di = 1.14 and Di = 6.99; = (2.0, -2.0, 0.0)' with 
d 2 = 0.224 and D 3 = 1.00; and x% = (2.1, -2.1, -0.2)' with d 3 = 0.0403 and D 3 = 0.14. This example 
can be found in |13j . 

4.2.6 Example 6 

A system of four equations F@ (x) — found in [3] , and given by 

yz + t(y + z) = 0, 
xz + t(x + z) = 0, 
xy + t(x + y) = 0, 
xy + xz + yz = 1, 

is solved. Two roots of F 6 (x) — are (x,y,z,t) = ± , , — ^— 1. Here, we study the 

\ 3 3 3 6 / 
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Table 4: Numerical results for system F±{x) = . 



















coc 










ECOC 










ACOC 








Method 


Xo 


k 




Xk) 


'IL 


IK 


-1 


Hoc 


Ap, 






II 






Ap fc 




II 






Apt 




NM 


Xq 


12 


8.06 • 


10" 


-4464 


3.46 ■ 


10 


-2232 


1.54 • 


10" 


-5 


4.92 


■ 10" 


3348 


1.50 


• 10" 


-3 


1.42 


■ 10" 


-1116 


1.54 ■ 


10" 


-5 


p = 2 


t 2 


12 


2.82 • 


10" 


-3616 


2.00 ■ 


10 


-1808 


6.03 • 


10" 


-5 


6.33 


• 10" 


2712 


9.55 ■ 


10- 4 


■(*) 


9.85 


• 10 


-905 


6.03 ■ 


10" 


-5 




x 3 




11 


6.26 • 


10" 


-3517 


9.42 ■ 


10 


-1759 


3.03 • 


10" 


-4 


1.28 


■ 10" 


2637 


9.10 


■ 10" 


-4 


7.31 


■ 10 


-880 


3.03 ■ 


10" 


-4 


AMN 




9 


1.58 • 


10" 


-3112 


7.37 ■ 


10 


-1038 


3.28 • 


10" 


-7 


1.11 


• 10" 


1728 


3.71 


• 10" 


-7 


2.48 


• 10 


-692 


3.28 ■ 


10" 


-7 


p = 3 


T 2 


8 


1.89 • 


10" 


-5422 


6.64 ■ 


10 


-1808 


9.64 • 


10" 


-5 


7.56 


■ 10" 


3013 


4.23 


■ 10" 


-4 


1.10 


■ 10" 


-1205 


9.64 ■ 


10" 


-5 




^0 


7 


1.30 • 


10" 


-3631 


6.64 ■ 


10 


-1211 


2.16 • 


10" 


-4 


3.59 


■ 10" 


2017 


2.51 ■ 


10" 5 


! W 


1.07 


■ 10 


-807 


2.16 ■ 


10" 


-4 


HMN 


^0 


9 


3.52 • 


10" 


-4442 


5.29 ■ 


10 


-1481 


5.23 • 


10" 


-4 


1.51 


■ 10" 


2466 


6.13 


■ 10" 


-3 


5.38 


• 10 


-988 


5.23 ■ 


10" 


-4 


p = 3 


x 2 


7 


8.20 • 


10" 


-7487 


8.96 ■ 


10 


-2496 


6.96 • 


10" 




1.63 


■ 10" 


4159 


1.33 


■ 10" 


-3 


1.82 


■ 10" 


-1664 


6.99 ■ 


10" 


-5 




Xq 


7 


1.67 • 


10" 


-4248 


1.90 ■ 


10 


-2830 


4.03 • 


10" 


-4 


2.16 


■ 10" 


2360 


4.37 


■ 10" 


-3 


6.83 


■ 10" 


-2830 


4.03 ■ 


10" 


-4 


FDN 


■'- o 


8 


2.19 • 


10" 


-8244 


1.60 ■ 


10 


-2748 


1.92 • 


10" 


-4 


1.96 


• 10" 


4580 


1.20 


• 10" 


-3 


1.23 


■ 10" 


-1832 


1.92 ■ 


10" 


-4 


p = 3 


X 2 




8 


9.95 • 


10" 


-4626 


2.08 ■ 


10 


-1542 


1.40 • 


10" 


-3 


3.35 


• 10" 


2570 


6.59 


• 10" 


-3 


1.61 


■ 10" 


-1028 


1.40 ■ 


10" 


-4 






7 


4.67 • 


10" 


-3419 


3.32 ■ 


10 


-1140 


1.39 • 


10" 


-3 


1.46 


■ 10" 


-633 


8.44 


• 10" 


-4 


1.05 


■ 10" 


-1899 


1.39 ■ 


10" 


-3 



convergence of the methods towards 

a = — « (0.5773502692, 0.5773502692, 0.5773502692, -0.2886751346)* . 

Initial conditions are = (0.5, 0.5, 0.5, 0.2)* with d x = 0.4887 and D x = 0.45; x 2 = (0.55, 0.55, 0.55, -0.1)* 
with d 2 = 0.1887 and D 2 = 0.1925; and x% = (0.6, 0.6, 0.6, -0.3)* with d 3 = 0.02265 and D 3 = 0.08. Like 
in the second example HMN method is of 4th order (*). In Table [H] in the methods marked (*) j = 5 was 
used. 

4.2.7 Example 7 

In [7] the following boundary value problem 

j/" + y 3 = 0, y(0) = 0, y(l) = l, 

was posed. To solve it, we consider the following partition of the interval [0, 1] : 

Uq = < iti < 112 < ■ ■ ■ < u n -i < u n = 1 , Uj + i = Uj + h, h = 1/n . 

Let us define y = y(u ) =0, yi — y{u\) , . . . ,y n ~i = y{u n -i) , y n = y{u n ) = 1 . If we discretize the 
problem by using the following numerical formula for the second derivative 

ri'» W - 1 ~ 2 g + 2W * 1 , k = 2,3,...,n-l, 

we obtain a (n — 1) x (n — 1) system of nonlinear equations: 

2yi - h 2 y\ - y 2 = , 
-y k -i + 2y fc - h 2 yl - y k +i = 0, k = 2, 3, . . . , n - 1 , 
y„_2 + 2y„_i - tfyl^ = 1 , 
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In particular, we solve this problem for n = 10. The solution a is 

(.105541119905921, .211070483662496, .316505813937525, .421624081569127, .525992841283953, 
.628906344657317, .729332377591977, .825878904047790, .916792309006097)* . 

We start with the points Xq,Xq,Xq, defined by = (1,0,-1,0,1,0,-1,0,1)* with di = 1.7290 
and £>i = 1.990; x 2 = (0, 0, 0, .5, .5, .5, 1, 1, 1)* with d 2 = 0.3165 and D 2 = 0.5012 and = 
(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9)' with d 3 = 0.02933 and D 3 = 0.007290. In Table we used j = 5 
in MN method (*). 

Finally, after performing the numerical computation of the previous examples, we have realized that there 
are not significant changes in the components of the approximations of the error vectors when using the 
above defined norm. 

5 Concluding remarks 

A generalization to several variables of a technique used to compute analytically the error equation of iterative 
methods without memory for one variable is presented. The key idea is to use formal power series. 

So far, using the definition of COC for one variable, an iterative method is applied to a set of functions and 
numerical verification of the value of local order claimed in the error equation is performed. In this paper we 
generalize this procedure to several variables. Furthermore, when the root is not known, as usually occurs in 
real problems, we have overcome this situation introducing the numerical computation of the order ECOC 
and ACOC. 

To obtain an approximation of the root with high precision it is necessary to define an adaptive float 
arithmetics allowing us, in each step, to get the appropriate length of the mantissa. For several variables 
this increasing value of the mantissa is defined according to the knowledge or not of the root. 

To illustrate the technique presented seven examples are worked and completely solved. In each one, four 
iterative methods have been carried out. Explicitly, agreements and differences are pointed out. 

After comparing COC, ECOC and ACOC we think that ACOC is the best when the root is known or not 
due to its speed of convergence and its adaptability to the definition of the mantissa length. 
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Table 5: Numerical results for system F§(x) = . 



coc 



ECOC 



ACOC 



Method xq 



11^(^)11 



Ap fc _x 



\\Sh\\ 



NM x 
P = 2 xl 

^3 



14 5.21 • 10" 



12 1.94-10" 



11 7.33-10" 



1.89 • 10" 1490 
2.17- 10- 1970 
6.97 • 10~ 1931 



1.01 • 10" 4 
7.69 • 10" 4 
1.57 • 10" 3 



1.81 ■ lO -2235 
5.85 • 10~ 2955 
5.28 • io -2896 



1.47 • 10~ 4 
1.23 • 10~ 3 
8.40 • 10~ 4 



9.60 • 10" 
5.21 • 10" 
7.48 • 10" 



1.01 ■ 10" 4 
7.69 ■ 10~ 4 
1.57 ■ 10~ 3 



AMN x\ 
P = 3 xl 
x 



9 6.58 • 10" 34S1 
8 5.97 - IO" 6440 
7 2.70 • 10" 4280 



1.02 ■ io- 1160 
4.99 • 10~ 2147 
2.98 • 10~ 1427 



1.63 • 10" 
6.63 • 10" 
1.67 • 10" 



2.75 ■ IO" 1934 
9.62 • 10" 3578 
2.27 ■ io" 2378 



3.42 • 10~ 4 
1.27 • 10~ 3 
2.74 • 10~ 3 



2.71 • IO" 774 
7.16 • IO" 1432 
4.33 • 10" 952 



1.63 ■ 10" 
6.63 ■ 10" 
1.67 ■ 10" 



HMN xl 

o 2 

P = 3 x 

x 3 



8 5.99-10" 



1.67 • 10" 



7 8.00-10" 



1.37 ■ 10~ 925 
8.94 • icr 1022 
6.99 • 10" 1738 



6.60 • 10"" 15 
4.27 • 10" 15 
2.49 • IO" 17 



6.33 ■ 10" 



3.11 ■ IO" 1 "" 
9.60 ■ 10" 2897 



9.51 • 10~ 13 
1.09 ■ 10~ 12 
1.15 • IO" 13 



4.63 • 10" 



3.49 • 10" 



5.28 ■ 10~ 31 



1.70 ■ 10" 



1.37 -IO" 1159 8.81 ■ IO -35 



FDN 



P - 3 x 



16 1.21 • 10" 



8 3.71-10" 



4.46 ■ IO- 1001 
6.10 • IO -1758 
9.88 • 10~ 1294 



3.54 • 10" 5 



1.28 • 10" 



4.26 ■ icr 1DO ° 
1.47- 10" 2929 
8.35 • io" 2156 



7.99 ■ 10" 
9.54 • 10" 
1.26 • 10" 



9.55 • 10"°°° 
1.85 • 10" 1172 
7.40 • IO" 863 



3.54 ■ 10" 



1.28 ■ 10" 



Table 6: Numerical results for system F§(x) = . 



















coc 










ECOC 










ACOC 








Method 


x 


k 


\\F & (x k ) 


'Hoc 


\\ e k 


-i 


Hoc 


Ap 


k-l 




II 


e*lloo 


Ap fc 




IIMoo 


Ap fc 




NM 


Xq 


12 


3.20 • 


10" 


-4679 


2.96 ■ 


IO - 


-2339 


5.14 ■ 


10" 


-4 


2.28 


10 -3509 


7.21 


• 10 


-4 


8.96 • 


IO" 1172 


5.14 


■ 10 


-4 


p = 2 


Xq 


12 


3.75 • 


10" 


-3304 


1.07 • 


10 


-1650 


7.28 ■ 


10" 


-4 


2.52 


■ IO" 2477 


4.86 


• 10" 


-4 


2.35 


• IO" 827 


7.28 


■ 10 


-4 


(*) 


T 3 


11 


3.60 • 


10" 


-3514 


1.10 • 


10 


-1757 


1.24 • 


10" 


■15 


3.37 


2Q-2636 


8.96 • 


10" 


18 


3.08 


• IO" 879 


4.72 ■ 


10" 


■33 


AMN 


Xq 


8 


4.38 • 


10" 


-7501 


3.15 • 


10 


-2498 


1.72 


10" 


-3 


1.67 


IO" 4166 


1.03 


• 10" 


-3 


4.46 • 


^q-1668 


1.72 


■ 10 


-3 


p = 3 


Xq 


7 


1.16 • 


10" 


-3528 


2.17 • 


10 


-1174 


3.65 


10" 


-3 


7.51 


. jq-1959 


2.19 


• 10" 


-3 


3.47 


• IO" 785 


3.65 


■ 10 


-3 


(*) 


r 3 


7 


1.13 • 


10" 


-3752 


1.63 ■ 


10 


-1251 


1.17 • 


10" 


13 


2.06 


2o -2085 


1.32 • 


10" 


21 


1.26 


• IO" 834 


4.71 ■ 


10" 


■35 


HMN 


Xq 


6 


1.37 • 


10" 


-4681 


3.30 • 


10 


-1168 


8.24 ■ 


10" 


-3 


4.07 


. j^q-2046 


4.71 


• 10" 


-3 


1.23 


• IO" 878 


8.24 


■ 10 


-3 


p = 4 


Xq 


6 


1.30 • 


10" 


-6611 


1.07 • 


10 


-1650 


5.83 ■ 


10" 


-3 


1.73 


-^q-2890 


3.34 


• 10" 


-3 


1.61 • 


10 -1240 


5.83 


■ 10 


-3 


w 


Xq* 


6 


3.24 • 


10" 


-7028 


1.10 ■ 


10 


-1757 


2.53 • 


10" 


20 


1.87 


10" 3075 


2.40 ■ 


10" 


25 


1.71 ■ 


IO" 1318 


9.09 • 


10" 


-46 


FDN 


Xq 


8 


8.73 • 


10" 


-6401 


1.50 ■ 


10 


-2131 


2.01 


10" 


-3 


2.48 


. j^g-3554 


1.21 


• 10" 


-3 


1.65 • 


IO" 1423 


2.01 


■ 10 


-3 


p = 3 


Xq 


7 


9.61 • 


10" 


-3188 


7.57 • 


10 


-1061 


4.04 ■ 


10" 


-3 


1.62 


■ 10" 1769 


2.43 


• 10" 


-3 


2.15 • 


io- 709 


4.04 


■ 10 


-3 




r 3 

Xq 


7 


1.43 • 


10" 


-3432 


6.51 ■ 


10 


-1145 


3.72 • 


10" 


16 


1.20 


■ IO" 1907 


4.73 ■ 


10" 


19 


1.85 ■ 


10~ 763 


1.47 ■ 


10" 


■35 
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Table 7: Numerical results for system F- 7 (x) = . 









coc 




ECOC 




ACOC 




Method xo 


k 


ll*V(*k)lloo 


IK-llloo 




Ffclloo 


Ap fe 


Halloo 


Ap fc 



NM 
p = 2 
(*) 



: 12 6.95 • 10~ 3714 
x n 



7.42 • IO" 3010 
4 11 1.39 -10 - 4759 



1.92 
1.98 
2.72 



■^q-1856 
10 -1504 
10 -2379 



8.19 • IO -11 
4.50 • 10~ 9 
1.23 • 1CT 9 



1.06 
1.11 
5.65 



. 10 -2784 
jq-2256 
■j^g-3569 



1.02 ■ 10~ 10 
1.81 • 10~ 9 
1.51 • 10~ 9 



5.52- 10- 929 
5.61 • 10~ 753 
2.08 • 10~ 1190 



8.19 ■ 10 -11 
4.51 ■ 10~ 9 
1.23 ■ 10~ 9 



AMN 
p = 3 



HMN 
p = 3 



FDN 

p = 3 



1 



Xq 

r 3 



8 3.81 • io~ 4298 
7 2.35 • 10~ 3175 
7 



1.90 • 10~ 5349 



2.22 

4.05 
8.17 



10 -1432 
1Q - 1783 



5.53 • 10"" 
1.55 • 10~ 5 
1.86 • 10~ 5 



2.24 
1.32 
1.98 



jq-2387 

■ io- 1763 
io- 2971 



1.91 • 10~ 6 
1.18 • 10~ 5 
2.57 • 10~ 6 



1.01 • 10~ 955 
3.26 ■ 10~ 706 
2.42 • io~ 1189 



5.53 ■ 10~ 6 
1.55 ■ 10~ 5 
1.86 ■ 10~ 5 



:J 8 1.23 -IO -3764 1.35 -IO -1254 2.68 - 10~ 6 5.68 ■ 10~ 2091 1.56 • 10~ 6 4.20 • 10~ 837 2.68 ■ IO -6 
■I 7 2.12 -IO -3201 7.52 • io~ 1067 3.78 • 10~ 5 4.61 ■ 10~ 1778 4.38 • 10~ 5 6.13 • 10~ 712 3.78 ■ 10~ 5 
7 3.45 -10~ 4877 1.90 ■ 10~ 1625 3.45 • 10~ 6 4.67 • 10~ 2709 1.79 • 10~ 6 2.45 • 10~ 1084 3.45 ■ 10~ 6 



xl 8 9.33 - io~ 6553 5.38 -IO -2184 1.20 • 10~ 7 6.10 ■ 10~ 3640 2.81 • 10~ 7 1.13 • 10~ 1456 
xl 8 3.54 ■ io~ 3092 1.81 • IO -1030 1.47 -IO -7 2.14 ■ 10~ 1717 1.71 • IO" 6 1.18 • 10~ 687 
xl 7 6.47 - io~ 476D 2.21 • 10" 1586 3.62 • 10~ 8 6.43 ■ IO -2644 2.69 • 10~ 8 2.91 • 10~ 1058 



1.20 ■ 10~ 7 
1.46 ■ 10~ 7 
3.62 ■ 10~ 8 
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